Suppression of thermally activated escape by heating 
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The problem of thermally activated escape over a potential barrier is solved by means of path 
integrals for one-dimensional reaction dynamics with very general time-dependences. For a suitably 
chosen, but still quite simple static potential landscape, the net escape rate may be substantially 
reduced by temporally increasing the temperature above its unperturbed, constant level. 
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Thermally activated escape over potential barriers is 
of relevance in a large variety of physical, chemical, and 
biological contexts In the most common case, a po- 
tential barrier A[/ much larger than the thermal energy 
kT yields an escape rate exponentially small in bJJ jkT . 
A first major generalization, of importance for concep- 
tual reasons as well as due to numerous applications, are 
periodically modulated potentials 0, S [1, 0, (sl i resulting 
in a renormalization of A[/ which depends in a very com- 
plicated manner on the details of the model. In all these 
cases the rate is thus a very rapidly increasing function of 
the temperature T . In our present work, the main focus 
is on the complementary case of a static potential and 
a time dependent temperature 'T]. In particular, we 
will demonstrate that in a suitably chosen, but still fairly 
simple and generic potential landscape, the escape rate 
of the unperturbed system at constant temperature may 
decrease upon temporally increasing the temperature. In 
view of the above mentioned results for constant temper- 
ature, this is a quite unexpected and counter-intuitive re- 
sult. Indeed, given that thermal noise is indispensable to 
escape, one would expect that an "extra dose" of noise 
should always enhance escape. Somewhat reminiscent 
previous findings always concern quite different types of 
systems: Dissipative quantum tunneling in the deep cold 
[8|, activationless electron transfer [j^], complex protein 
dynamics near the solvent glass transition models 
without a barrier against deterministic escape [lll|, or 
non-dynamical systems [T2| . 

As our main tool, we put forward a new path inte- 
gral approach, which unifies and extends several related 
approximations 0, i, i, d, i, [H, [H. Briefly, in differ- 
ent parameter regimes of the temporal modulations, the 
most relevant escape paths are of quite different char- 
acter. Therefore, each regime was so far treated sep- 
arately and the crossover omitted. Here, all potentially 
relevant paths are represented in terms of a suitable, one- 
dimensional parametrization and are kept till the final 
rate formula via an integral over all of them. 

We consider the ovcrdamped ID Langcvin equation 



?7(t)i(<) = -U\x(t), t) -f ^27j{t)kT{t) m (1) 

with time-dependent friction r](t) > 0, temperature 
T(t) > 0, and potential U{x,t). Dot and prime indi- 
cate temporal and spatial derivatives, k is Boltzmann's 
constant, and thermal fluctuations are modeled as usual 



[ij by (5-correlated Gaussian noise ^(t). For T — ^ 0, 
the deterministic dynamics is required to exhibit exactly 
one stable orbit (attractor) Xs{t) and one unstable or- 
bit (repeller) > Xs{t). Our main interest con- 
cerns the noise induced transitions of x{t) across Xu{t) 
for small but finite temperatures T(t), quantified by the 
rate T{t) := —n{t) at which the probability n{t) that 
x(t) < Xu(t) changes in time. 

To avoid unnecessary complications, we focus on ini- 
tial conditions x{to) = Xs(to), and we require the exis- 
tence of D := limf^oo Jq ^Kt) with h{t) := kT{t)/r]{t). 
Next, we divide ([T]) by rjlt) and employ transformed times 
i{t) :— dTb{T)/D, positions x{i) := x{t{i)), and forces 
F{x,i) := -DU'{x,t{t))/kT{t{i)), yielding, after drop- 
ping again the tildes. 



x{t) = F{x{t),t) +V2D ^{t) 



(2) 



In the general formalism, we will work with ([2]), while 
specific examples will refer to ([T]). The corresponding 
(back-)transformation of the rates r{t) = h{t)T{t{t)) / D 
readily follows from the obvious transformation of the 
probabilities n{t) = n{tit)). 

We first recall some basics, previously derived and dis- 
cussed in detail in Refs. [H']: For any given initial con- 
dition x(io) = X(), the probability density to find the 
stochastic process ([2|) at any "final" time t/ > <o at the 
position Xf can be represented as path- integral 



piXf,tf \ xo,to) 



(to) 



Vxit) e-^[-Wl/^ , (3) 



with action S[x{t)] := J^^ dt [x{t) - F{x{t),t)]'^ /4. Once 
this formal integral is evaluated, the rate follows as 

T{t) = -Ddpixu{t),t\xsita),ta)/dxu{t) . (4) 

For small D, the integral ([3]) is dominated by the path 
q{t) which minimizes the action S[x{t)] and thus satisfies 
the Euler-Lagrange equation 

m = ~p{t) F'{q{t),t) , p{t) :- q{t) - F{q(t), t) (5) 

with boundary conditions (/(to) = and q{tf) — Xf. 
Accounting for all paths x{t) "close" to q{t) by means of 
a functional saddle point approximation in Q yields 



p{xf,tf\xo,to) ^ [AirDQitf)] ^e 



'S[q(t)]/D 



(6) 
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where Q{t) satisfies Q{to) = 0, Q{to) = 1, and 



Q{t) = -[2Q{t)F'{q{t),t)] 



Qit)p{t)F"{q{t),t) . (7) 



Eqs. (HI), ^ yield an approximation to the rate 
r(t), which in principle becomes asymptotically exact as 
D ^ for any given t > to with a unique absolute mini- 
mum of S'[a;(t)], which generically is the case. But under 
many circumstances of foremost interest (e.g. relatively 
large t — to) even fairly small D are still far from this 
asymptotic regime, i.e. the saddle point approximation 
([6]) does not properly account for all relevant paths in 
([3]): Basically, a typical escape path x{t) spends almost 
all its time near Xs {t) , then crosses over into the vicinity 
of Xu{t), and remains there for the rest of its time. Any 
other behavior would yield a much larger action ^[^(t)] 
and thus is negligible in ([3]). However, rather different 
cross over "time windows" may still lead to almost equal 
S[x{t)], and a simple saddle point approximation is un- 
able to properly account for such quite remote regions in 
path space. In some cases, there may exists further local 
minima of 5'[a;(t)] and additional saddle point approxi- 
mations around each of them may save the case [3, @. 
The remaining problem is to keep track of all relevant 
minima and not to double count their neighborhoods if 
they get too close in path space. In other cases, e.g. for 
i- independent 77, T, and [/ in ([1]), there is a continuous 
"soft direction" in path space, invalidating plain saddle 
point methods altogether p^ . 

To overcome these problems we impose on top of the 
boundary conditions x{to) = Xs{to) and x{tf) = Xu(tf) 
the extra condition that x{t) arrives at some intermedi- 
ate point Xi at a given time U, and in the end integrate 
over all ti £ [to,tf] The pertinent formal relation, 

satisfied by the conditional probability density ([31), is 



p{Xf,tf\XQ,tQ) 



dU p{xf,tf\xi, ti) ^xlU\xo, to) (8) 



where 'i'x{t | xq, ^o) denotes the first passage time density 
across x, given x{to) = xq. For simplicity only, we as- 
sume from now on that Xi is located well in between Xs (t) 
and Xu{t) and is independent. Then, all non-negligible 
paths in ([3]) starting from x{ti) = Xi must immediately 
cross over to Xu{t) and thus admit for p{xf,tf\xi,ti) in 
([8]) a saddle point approximation ([6]) free of all the above 
mentioned problems. Focusing on xq = Xs{tQ) accord- 
ing to (HI), an analogous approximation ^ holds for 
p{xi,ti\xo,to), since all relevant paths in ([3]) now may 
leave the vicinity of Xs{t) only in the very end. By defi- 
nition, p(xi,ti\xo,to)dxi is the probability that x(t) from 
([2]) is encountered within [xi,Xi + dxi] at time ti, given 
x{to) ^ xo = Xs{to). Most such x{t) closely resemble 
the most probable path q(t) connecting q{tQ) — xq with 
q(ti) = Xi. On the other hand '^x:iiti\xo,tQ)dti is the 
probability that x{t) crosses Xi for the first time dur- 
ing i e [ti,ti + dti] . It seems reasonable to guess that 



most such x{t) once again closely resemble q{t). Hence, 
p{xi,ti\xo,to)dxi will essentially account for the same 
"events" as xiiti\xo,tQ)dti provided we relate the con- 
sidered intervals dxi and dti via dxi ~ q{ti)dti. Up to 
finite- Z3 corrections we thus obtain 



'^xi{U\xo,to) = q{ti) p{xi,t^\xo,to) 



(9) 



where q{t) satisfies q{to) — xq — Xs(to), q{ti) — Xi, 
and (O. More rigorously, our key relation ^ follows 
by adapting Ref . [la] to evaluate the derivative by x/ of 
([5|) in the limit Xf Xi. Details will be given elsewhere. 

The evaluation of the escape rate ([![) by means of ([5]), 
([H) , and ^ is the first main result of our present work. 
Similarly as in [3, d, 0| , closed analytical solutions of 
the concomitant differential equation ([5]), ^ are only 
possible for special F{x,t). To this end, we focus on 
piecewise parabolic potentials U{x,t) in ([IJ, correspond- 
ing to piecewise linear force fields in ^ of the form 



F{x<0,t) = Xs{t){x-y,{t)) + fit) 
F{x>0,t) = Xu{t)ix^yuit)) + fit) 



(10) 



with As(i)?/s(t) = A„(i)2/„(i) (continuity at x = 0). Fur- 
ther, the existence of stable and unstable orbits with 
Xsit) < Xi and a;„(t) > Xi is required, in particular 
Usit), Xs{t) < 0, yu{t), A„(i) > 0. For the natural choice 
Xi = 0, a straightforward but somewhat tedious calcula- 
tion 0, @| then yields for the rate ^ the result 
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$(i,r, to) 

Zit,T, to) 
YiT,to) 
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Zit,T, to) ^^^(^t,T,to)/D 
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4IuiT,t) 4/,(T,io) 
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(11) 
(12) 

(13) 
(14) 

(15) 
(16) 



We have verified that previous findings for time- 
periodic d, H, @ and time-independent systerns ([1]) [l^ 
are recovered as special cases. Those from [4^ are for- 
mally similar but contain quantities (called £ and s((/))) 
which are not explicitly available in general. 

As a first example we consider the dynamics ([Ij with 
constant friction 77 (t) = 1, a temperature pulse T(i) ac- 
cording to Fig. la, and a static, piecewise parabolic po- 
tential Uix), see Fig. lb. Already for the moderately 
small temperatures from Fig. Ic, the accuracy of the an- 
alytical approximation ([TT|) is quite good. We found that 
it quickly improves even further upon decreasing tem- 
peratures. After initial transients (omitted in Fig. Ic), 
Kramers rate is recovered until the temperature pulse 
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FIG. 1: (a) Temperature pulse with parameters To, AT, and 
T. (b) Piecewise parabolic potential with a well at Xs = 0, & 
barrier at Xu = 1, curvatures U"{xs) = 4, U"{xu) ~ —4, and 
barrier height AU — U{xu) — U{xs) ~ 1. (c) Time dependent 
escape rates for kTo — kAT = 0.1, = 1, and r — 0.05, 
T = 0.1, r = 0.2, T = 0.4, r = 0.8 (bottom up). Squares: 
Precise numerical solutions of ^ with seed x{to = —5) — Xs- 
Solid lines: Analytical approximation (|lip for the equivalent 
transformed dynamics ((2}, (|10p . Dashed lines: Kramers rates 
0] for T = To (bottom) and T = To + AT (top). 
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sets in at t = 0. Then, the rate rapidly increases and ap- 
proaches the Kramers rate corresponding to Tb+AT, pro- 
vided the pulse lasts sufficiently long. Finally, an analo- 
gous relaxation back to the original Kramers rate follows. 
Discontinuities of T{t) entail jumps of r{t). While the 
initial transients are well understood [ij], to the best of 
our knowledge no previously existing analytical approx- 
imation would be able to faithfully describe the "per- 
turbed and interfering transients" for largely arbitrary 
pulses and pulse-sequences. 

Next, we consider ([IJ with a periodically pulsating 
temperature T{t) and a piecewise parabolic potential 
U{x) exhibiting two barriers and two wells, see Fig. 
2. Transitions from to described by the rate 

^i^2{t), those from a:| to by r2^i(i), and those from 
x| towards a; = cx3 by T2^oo{t)- After suitable time- 
and space-transformations (cf. each rate can be 

approximated according to (jlip . Since they are small 
(transitions are rare), it is sufficient - as far as the 
populations ni{t) and n2{t) of the two wells are con- 
cerned - to consider their averages over one period T, 
denoted by ri^2, ^2^1, and r2^oo- Then, the popu- 
lations n{t) := {ni{t),n2{t)) are governed by the master 
equation [l| n{t) = —Mn(t) with matrix elements Mn = 
-M21 = fi^2, M12 = -f2^i, and A/22 = f2_-^i-|-f2^oo- 
The smallest eigenvalue of M is denoted by T and repre- 
sents the ultimate rate of escape towards 00 after initial 
relaxation processes, governed by the other eigenvalue of 
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FIG. 2: (a) Periodic temperature pulses with parameters To, 
AT, T, and T. (b) Piecewise parabolic potential with wells at 
xf — 0, X2 = 2.04, barriers at Xi = 1.41, X2 = 3.6, curvatures 
U"{xl) = 1, U"{x'l) = -10, f/"(a;i) = 10, U"{x^) = -1, and 
barrier heights Af/1^2 = 0.9, = 1, AC/2^00 = 1.1. (c) 

Dependence of the effective rate F across X2 on r and AT /To 
in units of the rate f = 9.11 ■ 10"" for T{t) = To. Shown are 
analytical results for T = 30, kTo = 0.05, rj{t) = 1, obtained 
as detailed in the main text. 



M, have died out. Since the two eigenvalues differ by a 
huge, Boltzmann-Arrhenius-type factor, the total proba- 
bility ni{t) + n2{t) that x{t) < X2 is expected and numer- 
ically observed to actually exhibit a practically perfect 
exponential decay e~'"' for all t > 0. The analytical re- 
sults for r are depicted in Fig. 2c. Their agreement with 
our numerical findings for the decay rate (not shown) is 
comparable to Fig. Ic. 

The most striking feature of Fig. 2c is a substantial 
reduction of the net escape rate T upon superimposing 
temperature pulses of suitable duration r and ampli- 
tude AT to the "unperturbed" temperature Tq. Roughly 
speaking, T{t) in Fig. lb approaches the instantaneous 
Kramers rate the quicker, the larger the curvatures in 



4 



Fig. lb are [l^. Since the curvatures relevant for f2_+i 
are larger than those for f i_2 and f 2^00 (see Fig. 2b), 
sufficiently small r mainly affect r2^i and thus lead to 
a reduction of the net decay rate F. We verified that al- 
ready a single temperature pulse (Fig. la) indeed yields 
an analogous reduction of escapes events. Fig. 2c fur- 
ther shows that the effect is overruled by competing sec- 
ondary effects when r and/or AT become too small. 
Finally, we have obtained very similar results also for 
U{xl) < U{x2), but from the viewpoint of equilibrium 
rates [l|, the case U{xl) > J7(a;|) shown in Fig. 2 seems 
even more surprising to us. 

Experimentally, potentials like in Fig. 2b are ubiq- 
uitous in the context of chemical reactions. E.g. in the 
modified case U{xl) < U{x2) these are reactions proceed- 
ing in two-steps via an intermediate (metastable state 
2:2). Temperature pulses could be generated, among oth- 
ers 0, by means of a flashing black body radiator. More 
realistic are short laser pulses ]X3,, whose basic effects 
(on the reacting molecules and their environment) may 
still be roughly modeled by a temperature pulse. On the 
other hand, we expect that our above mentioned main 



finding will be qualitatively robust against various mod- 
ifications of the pulsed perturbation, including more re- 
alistic models for tailored laser pulses. A further exper- 
imental playground are colloidal particles in a suitably 
designed potential landscape by exploiting light [18,] , di- 
electrophoretic [l9[, or magnetic [2^ forces. Temperature 
pulses in the form of acousto-mechanical white noise may 
be generated by means of piezo elements [21]. 

Since temperature pulses in ([T]) are basically equivalent 
to potential modulations in ([2]), we arrive at yet another 
quite astonishing conclusion: The escape rate for a t- 
independent temperature and a potential landscape as 
in Fig. 2b may decrease if the amplitude (multiplicative 
factor) of the potential is temporally reduced without any 
other change of its "shape" . This effect should be readily 
observable with colloidal systems [ll,[l^,[23| and possibly 
also with cold atoms in laser induced optical lattices Q 
or in complex reaction networks (23j . 
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